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DIFFUSIVE AND RADIATIVE EFFECTS ON VAPORIZATION 


TIMES OF DROPS IN FILM BOILING 
by Kenneth J. Baumeister and Glen J. Schoessow T 


ABSTRACT 

Diffusive and radiative effects are incorporated into an 
analysis for the vaporization time of drops in film boiling. The 
momentum, energy, and continuity equations are solved with some 
appropriate simplifications so as to obtain a simple closed form 
solution for the overall film boiling heat transfer coefficient. 
Next, a theoretical expression for the droplet vaporization time 
is developed and compared to the measured vaporization times of 
water droplets vaporizing into air, argon, nitrogen, and helium. 
The agreement between experiment and theory is good. Under the 
helium blanket, the diffusive evaporative component is signifi- 
cant in comparison to the film boiling component. 

INTRODUCTION 

Studies of the vaporization time of liquid droplets in film 
boiling are basic to many technological problems. For example, 
current research in droplet film boiling on hot surfaces is con- 
cerned with such diverse problems as fuel vaporization on the 
j 
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manifold of the spark- ignition engine or the combustion chamber 
of the diesel and gas turbine engine (ref. 1); cooling or quench- 
ing of heated surfaces (refs. 2 and 3); liquid droplet removal in 
mist sections of boilers (ref. 4); and in areas of safe fuel 
handling where accidental contact of fuel with a hot surface is a 
risk. The recent hazard study (ref. 5) of the problems associated 
with the accidental spillage of liquefied natural gas (LNG) on 
water also relates to the vaporization of large liquid drops in 
film boiling. References 1 to 15 contain a comprehensive current 
literature summary of the recent work in droplet film boiling as 
well as a discussion of the phenomenon. 

In some of the above mentioned applications, drops vaporize 
in a saturated atmosphere of their own vapor, while in other 
applications the drop is exposed to an atmosphere in which diffu- 
sion can occur. With the exceptions of references 6 and 8, the 
vaporization times have been measured under conditions where the 
drop is exposed to air or other inert gases. Under such condi- 
tions, diffusion could significantly affect the total vaporiza- 
tion time. The purpose of the present paper is to obtain a 
closed-form solution for the effect of diffusion, as well as ra- 
diation, on the total vaporization time of drops in film boiling. 

In analyzing the experimental data, two models of the boil- 
ing process have been considered. The first model, see fig- 
ure 1(a), has no provisions for diffusion from the top of the 
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drop. In this model, the vapor is assumed to flow upward from 
beneath the drop and enshroud the upper surface. Under these 
conditions, diffusion from the upper surface can be neglected. 
Most of the analyses to date use this assumption with the excep- 
tion of reference 10 by Gottfried and Bell. Bell (ref. 7) sug- 
gested the model shown in figure 1(b). In this model, the vapor 
from underneath the drop disperses radially outward away from 
the surface thereby leaving the upper surface exposed to an air 
or inert gas atmosphere. Under these conditions, molecular dif- 
fusion could occur. The fact that the droplet's upper regions 
are usually a few degrees below saturation temperature tends to 
support this model. Furthermore, Wachters (ref. 8) observed that 
a droplet in air evaporates faster than a droplet in a saturated 
atmosphere, thereby indicating that diffusive mass transfer 
exists. 

In general, most of the data available are total vaporiza- 
tion times for pure liquid drops resting quiescently on a flat 
plate at high temperatures, exposed to an ambient air. Under 
these conditions, the diffusion contribution is generally small, 
less than 10 percent, compared to film boiling contribution. 
Consequently, it is difficult to discern which model is valid. 
Schoessow and Baumeister (ref. 6), however, performed an experi- 
ment where diffusion was a dominant factor. They noted that 
helium (see table 1) in association with water has a diffusion 
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coefficient which is a factor of four greater than an air-water 
combination. Figure 2 shows their data for the water droplet 
vaporization time in air, argon, nitrogen, and helium atmos- 
pheres. From this data and a dimensional analysis, Schoessow 
and Baumeister concluded that the second model (fig. 1(b)) is 
valid. To further support this conclusion, they presented a 
novel flow visualization technique which outlined the various 
domains associated with the second model. 

Gottfried and Bell (ref. 14) included diffusion and radia- 
tion in an exact numerical solution for the two region model 
shown in figure 1(b). Their predicted vaporization times for 
small spherical drops in film boiling were in good agreement 
with their experimental results. On the other hand, Baumeister, 
Hamill, and Schoessow (ref. 11) neglected diffusion and obtained 
a single closed form solution based on model shown in figure 1(a). 
Their analysis applies to all drop sizes, small spheres to large 
liquid spills. Their predicted heat transfer coefficients and 
vaporization times were also in good agreement with the experi- 
mental results of reference 14. However, in light of recent sub- 
stantiation of the two-region model, the simplified theory of 
reference 11 can only be used when diffusion is negligible. 

The numerical results of Gottfried and Bell (ref. 14) and 
the theory of Baumeister et al. (ref. 11) compare favorably when 
diffusion is not a dominant mechanism. Unfortunately, the numer- 
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ical results of reference 14 cannot be compared to the theory to 
be presented in this paper because the numerical results of ref- 
erence 14 applies only to specific drop sizes with an air atmo- 
sphere . 

In the present paper, diffusive as well as radiative contri- 
butions are incorporated into an analysis for the heat transfer 
coefficient and vaporization time of drops in film boiling. The 
momentum and energy equations are solved with some appropriate 
simplifications to obtain a closed solution. A closed solution 
was desired for the convenience of the engineering designer. The 
theoretical results are compared to the measured droplet vapori- 
zation times of reference 6 for water droplets vaporizing into 
air, argon, nitrogen, and helium. 

The analysis will treat radiation and diffusion as small 
perturbations to the well known solution of the Leidenfrost 
problem. The analysis is limited to two-component systems in- 
volving a pure single- component liquid with small solubilities 
for the gas and a single -component gas of high purity (with ex- 
ception of air). 
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METHOD OF ANALYSIS 

The experimental measured, vaporization time of a discrete 
liquid drop in film boiling can be determined by a direct inte- 
gration of an energy balance on the drop: 

-Ap L ^ = h t (V)A b (V)£T (1) 

where the total heat transfer coefficient, h T , and the drop's 
bottom area A B are dependent on the volume V of liquid that 
exists at any time t. Heat is assumed to reach the drop only 
through lower area of the drop. In equation (l), the total heat- 
transfer rate to the drop, represented by the right side of equa- 
tion (l), is set equal to the mass evaporation rate times the 
latent heat of vaporization. 

On the right-hand side of equation (l), the mechanism for 
energy transfer to the drop is assumed ( see appendix A) to be 
conduction across the vapor film (in creeping laminar flow) and 
radiation to the bottom of the drop. On the left-hand side of 
equation (l), the mass evaporation is composed of diffusive 
vaporization from the top of the drop and film boiling from the 
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bottom. 

The problem, of course, is to first relate A,, and h„ to 
the properties of the liquid and vapor, droplet volume, plate 
temperature, and the environmental conditions surrounding the 
drop. After this has been accomplished, equation (1) can be di- 
rectly integrated to determine the vaporization time. First, the 
geometric parameters including Ag will be discussed. Secondly, 
an expression for the heat transfer coefficient will be presented. 
Finally, the integration of equation (1) will be discussed. 

Geometric Parameter A„ 

D 

In this analysis, the drop is represented by a flat disk of 
radius r Q and height SL as shown in figure 3. The uniform 
vapor gap 6, shown in figure 3, is assumed to simplify the math- 
ematics. For the cases where diffusive evaporation is small, 
this model gives realistic predictions of the heat transfer coef- 
ficient and vaporization time (ref. 11). 

The relationship of Ag, r Q , and i to droplet volume are 
given in table 2. The geometric parameters were determined in 
reference 11 from numerical solutions for the drop shape. From 
the numerical results, three analytical approximations were ob- 
tained for the various volume ranges listed in table 2. 
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Heat Transfer Coefficient 


An energy balance on the droplet can be written in the fol- 
lowing form, 

kA AT AD MPA 

W 1 — + "raA 41 ' -O-X 4 )*. + 


B 


RT r 
s o 


On the left hand side of the above equation, the total heat trans- 
ferred, h^Ag AT, is broken up into a conduction and a radiation 
term. The kA^ AT/6 term in the above equation assumes that 
only conduction occurs beneath the drop with the associated lin- 
ear temperature profile. Convection heat transfer causes the 
temperature profile to deviate from the assumed linear value. 
Convection will be accounted for later by introducing the modi- 
fied latent heat of vaporization. On the right hand side of the 
above equation, the first term represents heat removed from the 
drop's lower surface by evaporation, while the second term rep- 
resents heat loss by diffusion from the upper surface. For this 
form of diffusive energy, the partial pressure of the liquid in 
the vapor environment is assumed small and an equivalent spheri- 
cal shape has been assumed. 

Dividing both sides of the above equation by Ag and AT 

gives 

P v xw w AD ab MP s A d 

"t ~ 6 rad " AT RT r AT A 

so B 

The axial vapor velocity w(6) at the lower surface of the drop 
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can be found by a solution of the momentum and energy equations 
in the vapor flow field beneath the drop. Substituting, the ex- 
pression for w(6) into the above equation leads after some alge- 
braic manipulation to an expression for the total heat transfer 
coefficient h^* 

The detailed steps of the derivation for h^ are given in 
the appendix of this report. For ease of presentation, only the 
final expression for the total heat transfer coefficient is pre- 
sented here in the body of the report. 

The total heat transfer coefficient h T can be expressed 

as: 


h T = 


jl. - 0.25 N*/V* 1/4 - 0.75 h rad ,.. o 


/h 


((A-26) designates equation in appendix) where 

3\l/4 A 


N 


H / Pr' 


0.68 S 


c \GH 


(2)-(A26) 


(A25) 


where the symbols are defined in the nomenclature. 

In equation (2) , h Q is the expression for the film boiling 
heat transfer coefficient when the diffusive contribution N 


(eq. (A25)) and the radiative contribution h , are zero. In 

rad 


this case, 


h 0 " 


* \l/4/k 3 X*(p L - p v )p v g^ 4 


ATyL 


(3)-(A15) 


The dimensionless volume and area parameter will be defined in 
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equations (9) and (10) . 

For convenience in the integration of equation (1) , 
tion (3) is rearranged in the following form 



equa- 


( 4 ) 


where h is a dimensionless heat transfer coefficient. In a 
o 

similar manner, the expression for the total heat transfer coef- 
ficient can be nondimensionalized by introducing equation (4) in- 


to equation (2). The expression for h^, becomes 

* 


* 

h T =~ 


h 


1. - 0.25 


N 


*1/4 


- 0.75 


rad 


where 


* 

h T = /, 3, * 


Ac 3 >*(p l - P v )e T g\ 1/4 
V PTpL / 


£ (G PrH*) 1/4 
L r 


and 


rad 


rad /. 3 , * _\l/4 

/k X (p L - P v )P v g\ 


ATpL 


(5) 


( 6 ) 


(7) 


In the next section, we will substitute the relationships 
for the area A and the total heat transfer coefficient h T 
into equation (1) and integrate. This will lead to expressions 
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for the total vaporization time which can be compared to the ex- 
perimental data in figure 2. 


Analysis of Vaporization Times 


Defining 



and h T as given in equation (4) , equation (1) can be conveni- 
ently written in dimensionless form as: 

* * * * 

-dV = h T A dt (12) 

* 

In this case, A represents the dimensionless bottom area of the 
flat disk model shown in figure 3. Integrating equation (12) 
gives : 



(13) 


The integration of equation (13) for 


* 

V 


greater than 0.8 
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requires that the integration be broken up into two ranges since 

* 

the form of and A are different in each volume range. 

* 

Consequently, for V greater than 0.8, the integration is from 
* 

V to 0.8 and from 0.8 to zero. A similar approach is used for 
& 

V greater than 155. 



(14) 


(15) 


I 


dV 


* * 

'155 h T A 


(16) 


where the appropriate form of the geometric parameter A and 
hj, is used in the evaluation of each of the above integrals. 

The above integrals are easily evaluated because all the integrals 
reduce to simple polynomials. 

The result of these integrations is displayed in table 3. 

a 

The formula relates the dimensionless time t to a pseudo liquid 
volume V + . When diffusion and radiation are zero the pseudo 


volume becomes equal to the actual volume. 
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COMPARISON OF EXPERIMENT TO THEORY 


The experiment data of reference 6 are compared to the 
theoretical predictions in figure 4. In evaluating N , the 
diffusion coefficients in table 1 are used, the fluid and vapor 
properties are evaluated at the film temperature and the ratio 
of diffusive to film boiling area (A,/A_) in N is taken as 
1.35. Some diffusive evaporation has been assumed to take place 
from the sides of the drop, hence the thirty-five percent increase 
in the diffusive area as compared to the film boiling area. This 
average area ratio applies to drops in the dimensionless volume 
range of 0.8 to 155, which is the range of the experimental data. 

Theory and experiment are in good agreement as seen in fig- 
ure 4. The data fall within a standard deviation of 5 percent. 

In comparison, if the diffusive component of evaporation is 
neglected (D^g = 0) , the helium data fall approximately 60 per- 
cent below the theoretical curve. 


LIMITATION OF THEORY 

Various simplifications (linearizations) were used in deriv- 
ing the expressions for the total heat transfer coefficient, see 
equations (A19) , (A21) , (A22), and (A23) . These linearizations 
lead to a singularity in equation (2) in the limit when 


0.25 N 
„*l/4 


rad 


+ 1 


+ 0.75 
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that is, the heat transfer coefficient becomes infinite. 

For water drops vaporizing into a helium atmosphere, this 
singularity occurs at a water droplet volume of approximately 
0.0001 cc. For large liquid drops, this singularity will have 
a negligible effect on the predicted total vaporization time; 
however, this will not be the case for very small liquid drop- 
lets. Consequently, it is suggested that the correlations in 
table III be used only in the range tested by experiment. From 
the present results, the last two equations in table 3 
(V >0.8) can be safely used. 

* 

For very small liquid volumes (V < 0.8), the numerical pro- 
cedure of Gottfried, Lee, and Bell in reference 10 should be used 
for predicting the droplet vaporization time. 

CONCLUSIONS 

A closed-form solution for the heat transfer coefficient and 
droplet vaporization time is in good agreement with data in which 
diffusion plays an important role. 

NOMENCLATURE 

bottom area of drop 

area of drop in which diffusion occurs 
dimensionless area, defined by eq. (10) 
specific heat at constant pressure of vapor 


B 


A 

C 
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D diffusion coefficient (see table I) 

A B 

3 2 

G defined by p (p_ - p )gL /y 

g coefficient of gravity 

g^ gravitational constant 

H defined by X/C AT 

P 

H defined by X /C^ AT 

h heat transfer coefficient 

film boiling heat transfer coefficient when radiation and 
diffusion are assumed zero 

— 4 4 / 

h , radiative heat transfer coefficient e T o(T - T )/(T - T ) 

rad L w s / w s 

total heat transfer coefficient to bottom of the drop 

k thermal conductivity of vapor 

L characteristic length, equation (11) 

& thickness of flat disk drop model 

M molecular weight 

N D £ dimensionless diffusion parameter defined by equation (A24) 

* 

N dimensionless diffusion parameter defined by equation (A25) 

P pressure 

P q ambient pressure 

Pr Prandtl number C y/k 

P 

P g saturation pressure 

R gas constant 

r radial coordinate 

r^ effective droplet radius 
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£ 

r dimensionless radius, r /L 

o o 

5 modified Schmidt number, y/D _(MP /RT ) 

C S S 

T temperature 

film temperature, equation (Al) 

Tp plate temperature 

T g saturation temperature 

AT temperature difference, (T - T g ) 

t time 

Jc 

t dimensionless time, equation (8) 

u radial velocity 

V droplet volume 

V dimensionless droplet volume equation (9) 

V + pseudo volume, see table 3 

w axial velocity 

w(6) axial velocity at bottom of drop 
Z axial coordinate 

6 vapor gap thickness 

vapor gap thickness when diffusion and radiation are not 
present 

emissivity of liquid 
latent heat of vaporization 


X 


17 


A* 


M- 

PL 

P v 

a 

a 


modified latent heat of vaporization (ref. 13) 


A 


< 2 


A 


> 2 


/ 7 Cp ar 

V 3 

\ 20 A ^ 


0.874 In ^1 + 


2A J 

C p & 


2A 




viscosity of vapor 
liquid density 
vapor density 
surface tension 
Stefan-Boltzmann constant 
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APPENDIX - HEAT TRANSFER COEFFICIENT 


The heat transfer coefficient from the plate to the drop is 
obtained by solving the momentum and energy equations for flow 
and energy transfer in the vapor gap beneath the drop. The 
physical properties are assumed constant and evaluated at the 
film temperature: 


T + T 



For creeping laminar flow with only heat conduction across the 
vapor gap, the governing differential equations are as follows: 
Momentum : 


0 = 


~ g c 9P. 
P v 9r 




(A2) 


0 = 


-s. 


ap 

9z 




Continuity: 


8u , u , 3w 

i7 + 7 + 3l 


= o 


(A3) 


(A4) 


Energy: 

h = | (A5) 

Boundary conditions: 


z = 0 

u = 0 

u = 0 

T = T 

P 

(A6) 

Z = 5 

u = 0 

w = w(6) 

T _ T 

s 

(A7) 

r = 0 

u = 0 



(A8) 
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(A9) 


(A10) 


(All) 


The set of equations presented above are similar to those 
given in reference 12 except that diffusion and radiation terms 
are added to the droplet energy balance. On the right-hand side 
of equation (All), the first term represents heat removed from 
the drop by evaporation from the lower surface, while the second 
term represents heat loss by diffusion from the upper surface. 
For the form of diffusive energy loss, the partial pressure of 
the liquid in the vapor environment is assumed small and an 
equivalent spherical shape has been assumed. 

The momentum equation along with the appropriate boundary 


conditions was solved in reference 12 to give 


... 2 , (p L - 

" (s) ' 3 Jl 


where V and A are dimensionless volume and area of the 


bottom of the flat disk model, respectively. They are defined 
by equations (9) and (10) in the body of this report. 
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Equation (A12) could now be combined with the interface 
energy balance, equation (All), to determine the gap thickness 6. 
Next, the heat transfer coefficient can be found by combining the 
gap thickness <5 with equation (A5) . However, for convenience, 
the vapor gap thickness and heat transfer coefficient will first 
be found for the case of pure film boiling. 

For pure film boiling, no radiation or diffusing, the inter- 
face energy balance becomes 


V " " P v Xw(6 o ) CA13) 

o 

where the subscript "0" has been added to signify that this is 
the gap thickness in the absence of diffusion or radiation. Sub- 
stituting the expression for w(S) from equation (A12) into equa- 
tion (A13) and solving for the vapor gap thickness gives 


= ( 3k atul y/y/iy 74 

o \2ttA(p, - P Jog] 


(A14) 


J L K v /v W W / 

From the definition of the heat transfer coefficient in equa- 


tion (A5) , 


h 

o 



1.21 


ATyL 


P v )p v g 



(A15) 


where the latent heat of vaporization X has been replaced by 

* 

the modified latent heat of vaporization X . The modified 
latent heat of vaporization accounts for the convective terms in 
the energy equation which were neglected in equation (A5) . The 
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function form of X was derived in references 12 and 13 and is 
given in the list of symbols. 

Now, substituting equations (A12) and (A13) into equation 
(All) gives 


X i P D AB Mf S / A d\ 

'C ' k l R Vo iT V A B / ' 


rad 


(A16) 


or multiplying the diffusive and radiative terms by the identity 


6 0 /$ 0 gives 


6 4 = 6 4 
o 


1 - N. 


rad \ 6 


DC 


h ; 5 
o / o J 


(A17) 


where 


„ A d H (pr 3 V /Z| 1 

a rtm 1 Am S \ ^ J 


DC A„ RT r h AT 
B s o o 


(A18) 


J c \G H / h*r* A B 

r o o 

* * 

The definitions of h and r are given in table 2. This 

o o 

dimensionless group is a measure of the importance of diffusion 
from the top of the drop to conduction below. 

Equation (A16) is the biquadratic equation for which solu- 
tions exist (ref. 16). Unfortunately, the more general solution 
does not lead to a closed-form solution for the total vaporiza- 
tion time. However, for the range of and h^^ associated 

with typical experiments, a simple iterative solution can be found 
by assuming 6 x 6^ on the right-hand side of equation (A17). 


This gives 


6 = yi - (n dc - h rad /h o )] 


i/a 


(A19) 
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The total heat transfer coefficient to the drop can now be 
found by substituting 5 into the left-hand side of equa- 
tion (All) . 

h 


t 1 - < N DC - h rad /h o>] 


1/4 + ^rad 


(A20) 


This form of the total heat transfer, although relatively simple, 
will still not lead to a closed-form solution for the total va- 
porization time. Therefore, equation (A20) is now written in the 
form of a binomial expansion in which only the first order terms 
are kept. The details are as follows. 

h 


T 1 - I ®DC - WV 


+ h 


rad 


(A21) 


or 


h T = h o 


1 + 4 N DC 


rad 


+ h j = h 
rad o 


3h 


1 + 


n dc + h 


rad 


(A22) 


Finally, inverting gives 

h 


T 1 - j- (N n _ + 3h ,/h ) 
4 DC rad o 


(A23) 


The diffusion parameter N D( , is a function of the dimen- 

A A 

s ionless volume V which is embedded in the terms h and r 

o o 

£ ^ 

in equation (A18) . Substituting in the value of h 0 r 0 from 
table 2 into equation (A17) gives 


N. 



DC 0.68 S 


(A24) 
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Defining 


N 


_1 H_ [ Pr 

°' 68 S c \G H* 
>• r 


3 \l/4 


\hj 


(A25) 


The heat transfer coefficient to the drop becomes 


1. - 0.25 N*/V* 1/4 - 0.75 h rad /h Q J 


(A26) 


In this case, N is independent of drop volume. This form of 
the total heat transfer coefficient leads to a closed solution 
for the total vaporization time. 
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TABLE 1. - DIFFUSION COEFFICIENTS (373 K 


AND ATMOSPHERIC PRESSURE) 


Liquid-gas 

combination 

d ab> 

cm^/sec 

Refer- 

ence 

Comments 

Water -air 

0.353 

16 

— 

Water-argon 

0.41 

17 

(a) 

Water-helium 

1.22 

17 

(b) 

Water-nitrogen 

0.40 

17 

(b) 


calculated by Slattery-Bird equation cor- 
rected for water (eq. (11-27)) as given in 
ref. 17. 


^Experimental data given in table 11-4 of ref. 17 
was linearly extrapolated to 373 K. 










TABLE 2. - RADIUS AND CONDUCTIVE HEAT TRANSFER COEFFICIENT FOR DROPLET FILM BOILING 




















TABLE 3. - THEORETICAL EXPRESSIONS FOR DROPLET VAPORIZATION TIME 


Dimensionless 

volume, 

* V 

v "7 

Drop shape 

v + 

* 

t 

/ 

V* < 0.8 

Small spheroid 

• 

V + - ^v* 5/12 - 0.624 N*V* 1/6 - 0.142 h* ad V* 1/2 ^ 

t* - 1.205 v +5/12 

VWWWWW^ 

0.8 < V* < 155. 

Large drop 

V + - ^v* 1/3 - 1.345 N*v* 1/12 + 0.995 N* - 0.465 h* ad V* 1/2 + 0.345 h*^ 3 

t* - 2.23 V +1/3 - 0.97 


V* > 155. 

Extended drop 
^Constant \ 

^ thickness J 

v + - | V* 1/4 + 0.0625 N* + 2.29 h* ad V* 1/2 - 0.524 N* - 32.3 h*J 4 

t* - 4.52 V +1/4 - 5. 

lt\\\\\\Y\\\YUH 


’.-.V -• 

-y ilNERT GAS . • 


SATURATED VAPOR 




•INERT GAS* 


FREE CONVECTION -DIFFUSION 
DOMINATED REGION 

v X/ ^ J J 

. ^ REGION 2^ 


't o-; 



REGION 1 

SATURATED 

VAPOR 


(b) DISPERSED ATMOSPHERE. 


Figure L - Velocity flow models for drops in film boiling. 
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Figure 2. - Total vaporization data for various water and 
inert gas combinations. 



Figure 3. - Schematic model for evaporation of flat disk. 


Figure 4. - Total vaporization time correlation for water 
droplets (data from ref. 6). 


NASA'L«wis-Com'l 



